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SUMMARY 

Various weighting methods are applied to typical nonuniform duct 
flow profiles to determine average flow properties. The analysis covers 
a range of subsonic duct Mach nunibers, but is confined to flows having 
uniform static pressure and total temperature. 

An averaging method is developed which yields uniform properties 
that reproduce the mass and momentum of the nonuniform flow. In con- 
trast, it is shown that the use of conventional weighting methods may 
result in large errors in these properties. These errors are shown to 
have varying significance depending on the applications to which the data 
are applied. 

It is also shown that nonuniform flows through variable-area duct 
passages will cause changes in average flow properties that are not as- 
sociated with the real thermodynamic flow path. 


INTRODUCTION 

In most calculations involving duct air -flow properties, it is not 
convenient to consider local flow variations within the duct. Therefore, 
the properties of the flow are treated as though they were uniformly 
distributed, and one-dimensional equations are applied to this uniform 
flow. Inasmuch as the real flow seldom approaches uniformity at planes 
of interest , the equivalent uniform flow must be determined by some 
method of averaging the properties of the real flow. 

This report presents the results of an analytical study made to de- 
termine the accuracy with which several commonly used averaging or 
weighting methods reproduce the real flow properties. The significance 
of inherent errors is illustrated for several common applications of 
duct flow data. Errors introduced through the application of one- 
dimensional relations to the uniform flow are briefly examined. 
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The study considers several typical velocity gradients but is con- 
fined to subsonic compressible flows with uniform static pressures and 
stagnation temperatures . 

(Since the present analysis was completed^ it has been found that a 
more generalized, qualitative analysis of the same problem is contained 
in ref. 1. ) 


ANALYSIS 

A uniform flow representing the flow properties of a nonuniform 
duct flow should satisfy the total energy, mass, and momentum of the 
real flow. For the special case considered herein in which the flow is 
assumed to arise from a uniform temperature source and to flow adia- 
batically to the measuring station, the total energy of the real flow 
can be reproduced by the assumption of constant total temperat*ure in 
the uniform flow at the source value. The determination of a uniform 
flow that will simultaneously satisfy the mass flow and the momentum in 
the real flow is more difficult. 


Mass -Momentum Method 


For the special case in which the static pressure and total 
perature are constant across the duct, the mass flow is given by 
equation 


m 



r\ 
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tem- 
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( 1 ) 


where M is the axial component of the local duct Mach number. (All 
symbols are defined in appendix A. ) 


In order for the mass flow in the representation to equal this in- 
tegrated mass flow, the uniform flow must satisfy the relation 
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( 2 ) 


where pg and M^ are the effective static pressure and Mach number, 
respectively. 


The integrated momentum of the real flow can be expressed by 

cp = p J (l + rM^)dA 


( 3 ) 
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Thus, the effective static pressure and Mach number must also satisfy 
the relation 


<P = Pe(l + tm|)A (4) 

By combining equations (l) to (4), the expression for the effective 
Mach number required to satisfy the total energy, mass flow, and momen- 
tum of the real flow becomes 


^ 

9 1 + 

where m and <p are integrated values determined from equations (l) 
and ( 3 ). 

Although for this analysis the static pressure is assumed to be 
constant across the real duct flow, this measured value of pressure 
cannot be used in conjunction with the effective Mach number determined 
from equation (S) to satisfy the real flow properties. Instead, a new 
effective static pressure must be determined from either the momentum 
or the mass -flow equations as 



( 6 ) 


This effective static pressure is never identical to the measured pres- 
sure if velocity gradients are present in the real flow. 

To complete the definition of the equivalent uniform flow, an ef- 
fective total pressure can be determined from the expression 

The flow quantities defined hy this method of averaging would he 
those obtained by mixing the measured profile to a uniform flow in a 
constant-area duct without wall friction. Mixing losses are inherently 
contained in the average flow quantities. 
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Conventional Weighting Methods 

The weighting or averaging methods commonly used to obtain uniform 
flow representations of nonuniform duct flows require either less com- 
plicated data-collection methods or less tedious calcxilation techniques 
than does the exact weighting procedure. Such methods result in in- 
herent errors in the representation of one or more of the basic proper- 
ties of the real flow. The required assumptions and applicable equa- 
tions for three of the more commonly utilized methods fol^l-ow. 

Mass-derived method . - When the mass flow in the duct is known from 
some independent measvirement, the measured static pressure at a station 
can be used in conjunction with the geometrical flow area A to define 
a uniform duct Mach niamber M^ that satisfies the mass flow by the 

equation 


M. 




M^ 
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2 “ 
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pA 


( 8 ) 


From this average Mach ntimber and the measured static pressxire , an 
average total pressure can be calculated as 


Pc = p(i + ^ 4 ) 


L- 

r-i 


(9) 


The momentm calculated from the measvired static pressure and the 
average Mach number becomes 


<p^ = p(l + tm2)a 


( 10 ) 


It is evident that the mass flow and total energy of the real duct 
flow are inherently satisfied by the mass-derived method of determining 
an average flow. There is no attempt in this method^ however, to satis- 
fy the momentum of the real flow. 


Mass -flow -weighting method , - A pitot -static survey of the flow at 
the desired duct station is frequently enployed to determine an average 
uniform flow. If it is assumed that the measured nonuniform flow can 
be brought to rest without mixing losses, the resultant pressure can be 
determined from the equation 


Pc = 




Ppv dA 


dm j’^'pV dA 


(lla) 
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For the special case in which the static pressure and total ten 5 )erat\xre 
are constant across the duct, the compressible form of equation (lla) 
becomes 
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(11b) 


The mass flow and momentum of the uniform flow having a total pres- 
sure defined by equations (ll) are not unique values. Their magnitudes 
depend upon the nature of additional assvimptions about the properties 
of the uniform flow. 


The measured static pressure at the duct station is often assumed 
to be the static pressure of the average flow. With this ass\m 5 >tion, a 
uniform duct Mach number can be defined by the relation 


1 



's 

The momentum for this uniform flow is given by equiation (lO). 
calculated mass flow becomes 


( 12 ) 


The 


m„ 



(13) 


The mass flow determined from equation (l3) will not correspond to the 
integrated mass flow which was used to determine the average total pres- 
sure in equation (lla). 


This anomaly between the integrated and calculated mass flows can 
be avoided by defining an average static pressuire p,, which, when used 

with the average total pressure from eqioations (ll) , will satisfy the 
integrated mass flow. The average Mach mmiber required to satisfy the 
mass flow under these conditions is given by 


M. 
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( 14 ) 
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and the resultant static pressure becomes 


Pc 


P 


c 



r 


r-i 

2 



(15) 


The momentum calculated from equation (lO) with either the measured 
or calculated v'alues of static pressure and corresponding Mach n'umber 
will not equal the integrated momentm. Equations similar to (l4) and 
(is) can he determined which would yield a static pressure and Mach num- 
ber for the uniform flow that would satisfy the real flow momentum. 

These flow properties would not satisfy the mass flow, however, and are 
not conventionally eir 5 )loyed. 


Area-weighting method . - When pitot-static flow surveys are em- 
ployed, the con 5 )lications of the calculation procedure can be reduced 
by using an area-weighted average total pressure determined from the 
equation 



(16) 


The remaining properties of the uniform flow are calculated by the 
eqviations used with the mass-flow -weighting method. As in the former 
method, several solutions for these properties are possible. Generally, 
the static pressure is assumed equal to the measured value. If inde- 
pendent mass measurements are available, a static pressure may be cal- 
culated to satisfy the mass flow. With compressible duct flow, the in- 
tegrated momentum will not be satisfied with either assumption. (For 
the incompressible case, a uniform flow defined by the total pressure 
from eq. (l6) and the measured static pressure will duplicate the real 
flow momentum. ) 


NUMERICAL CALCULATIONS 

The uniform flow properties of three arbitrary duct profiles were 
calculated by the mass-momentum weighting procedure and by the conven- 
tional weighting methods discussed in ANALYSIS. For simplicity, the 
ducts were assumed square with symmetrical two-dimensional profiles. 
The profiles considered were : 

(a) A power profile described by 

1 

M = 


(17) 
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(b) A discontinuous. 


(c) A linear profile 


separation profile represented by 


0 < X < 0.1 

M = 0 

(18a) 

0.1 < X < 1.0 

M = K 

(18b) 

of the form 



M = K(0.2x + 

0.8) 

(19) 


Each profile vas evaluated for a range of values of K (correspond- 
ing to the maximum Mach nuniber at the duct centerline) from 0 to 1.0. 


Mass -momentum method . - Equation (l) was integrated for each of the 
profiles to determine the mass flow actually contained in the duct. The 
integrals for the power and linear profiles were approximated by series 
expansion. The resultant expressions for the mass flow (valid for 
K < l.O) were 



pA 


0.875K 

0.9K^1 
0.9K + 


+ 0.070K^ 


1 



0.0738K^ 


0.00292K^ + 0.00025K'^ 

(Separation profile) 

0. 00307 4K^ + 0.00026K’^ 


(Power profile) 
(20a) 


(20b) 


(Linear profile) 
(20c) 


The actual momentum with the assumed profiles was obtained by inte- 
grating equation (3) with the resultant expressions 


■^ = 1 + 1.08889K2 
pA 

(Power profile) 

(21a) 

= 1 + 1.26K^ 

(Separation profile) 

(21b) 

= 1 + 1.1387k2 

(Linear profile) 

(21c) 


Effective values of duct Mach mmiber, static pressure, and total 
pressxire were determined from equations (s), (S), and (7), respectively. 

Mass -derived method. - By using the values of mass flow from equa- 
tions ( 20 )^ the properties of the uniform flow were determined from 
equations (8) to (lO). 
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Mas s-flov -weighting method . - The product of the average total pres- 
sure and the mass flow was obtained from equation (lib) and modified to 
the form 

integrated, yielded the following expressions; 

0.56K^ + 0.14K^ + O.OieK"^ + O.OOOTK^ (Power profile) 

(23a) 

2\4 

—■] (Separation profile) (23b) 

= 0.9K + 0.5904K^ +0.1476K^ +0.01665K'^ + 0.00071K^ (Linear profile) 

(23c) 

The values of integrated mass flow from equations (20) were then used to 
obtain the average total pressure. Equations (l2) to (l5) were \ised, as 
appropriate, to determine the calciilated average properties of the flow. 

Area-weighting method . - The average total pressure was obtained 
from equation (16), which becomes 


Equation (22), when 




P ^ i 

— — 4"^ = 0.875K + 
P pA 


= 0.9K^ + 



The resultant expressions , 


after integration, were 


(24) 


h. 

P 


1 +0.5444K^ +0.1H4K^ +0.00942K® +0.000204K® - . . 
7 

/ 

0.1 + 0 . 9(1 + (Separation profile) 

1 +0.5693k2 +0.1177K^ +0.00988K^ + 0.0002 IK® - . . . 


(Power profile) 

(25a) 

(25b) 

(Linear profile) 

(25c) 
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The average properties of. the flow were calculated from equations 
(12) to (15). 


RESULTS 

The profiles assumed for the numerical analysis were chosen to 
represent conservative flow nonuniformities as compared with those often 
experimentally observed. Typical profiles are presented in figure 1. 
Each profile in this figure corresponds to an effective duct Mach number 
of 0.2 as determined by the mass -momentum method. In addition to the 
Mach number profiles, the accompanying total. -pressinre variations (for a 
constant duct static pressure) are presented in the form of the local 
incremental deviation in total pressure from the mean effective value 
determined by the mass -momentum method. 

At the duct centerline (x = l.O), the maximum total -pressure devia- 
tion occurred with the power profile. In this case the local total 
pressure exceeded the effective average value by about 1 percent. With 
the wall pressure (at x = O) used as an indication of the other ex- 
treme in total-pressure deviation, the separated profile gave a maximum 

deviation of less than 3^ percent below the effective value. For p\n:- 

poses of qiialitative comparison, the wall static pressure that would be 
observed for a uniform duct Mach nimober of 0.2 is indicated in the fig- 
ure. It can be concluded, therefore, that all the assumed profiles 
represent moderate flow distortions. As a consequence, the errors that 
will be shown to accoii5)any the various weighting techniques are less 
than might be expected for practical flow problems. 

Figxcre 2 compares the static pressures that wovild be measured for 
each of the assximed profiles with the corresponding effective static 
pressures determined by the mass -momentum method. It is seen that the 
measured static pressure in a duct having nonuniform velocities will 
always be less than the effective static pressure required to describe 
the integrated flow properties in the duct. 

The deviation between meas\ired and effective static pressures in- 
creases as the effective duct Mach nuinber increases for the assiimed 
profiles. This results from the inherent nature of the profile assun^)- 
tions, wherein the magnitude of the total -pressure variation across the 
duct increases as the maximum duct Mach nuinber K, and hence the effec- 
tive Mach number, increases. The curves terminate at a value of K = 1 
for each profile. It is interesting to note that an effective duct Mach 
number of unity, as defined by the mass -momentum method, cannot be 
achieved with any nonuniform duct flow, regardless of the value of K. 
This restriction arises from the fact that the mass flow with uniform 
sonic velocity is greater than the mass flow in a nonuniform stream of 
the same area, whether subsonic or supersonic. 
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The average total pressure determined by each of the weighting 
methods for the assumed profiles is compared with the effective vsilue 
of total pressure from the mass-momentum method in figure 3. These to- 
tal pressures are independent of any further assTjmptions regarding the 
average static pressure in the duct. 

Mixing losses are inherently included in determining the effective 
total pressure by the mass -momentum method. Since the mass -flow- 
weighting method assumed no mixing losses, this method always yields an 
average total pressure that is greater than the effective value. On the 
other hand, the mass-derived and area-weighting methods yield average 
total pressures that are lower than the effective value. The errors 
with all methods of averaging increase as the severity of the profile 
increases whether through an increase in the value of K, and hence of 
M^, or from the nature of the basic profile shape. 

For weighting methods in which the uniform Mach ninnber in the duct 
is determined from the calcxilated total pressure and the measured static 
pressure, the combined errors in static pressure (fig. 2) and calculated 
total pressure (fig. 3) might be expected to produce significant errors 
in Mach number. This expectation is confirmed by the cvirves of figure 
4, which show that all the weighting methods yielded calculated Mach 
nimbers that were greater than the corresponding effective Mach numbers 
determined by the mass -moment \jm method. 

Any errors in the determination of static pressure and Mach number 
for the uniform flow will reflect as errors in the calctilated mass flow 
and momentum. The magnitude of these errors is illustrated in figures 
5 and 6 for the profiles and weighting methods considered. 

Inasmuch as all the averaging methods used the measured static 
pressure in the calculation of mass flow in figure 5 and of momentum in 
figure 6, the ratio of measured to effective static pressure was iden- 
tical. The differences in the calculated values therefore arise from 
the differences in Mach number computed by the various methods. With 
the mass-derived method, the calculated Mach number exactly satisfied 
the measured mass flow when used with the measured static pressure. 

The l>4ach number ratios indicated in figure 4 for the mass-derived meth- 
od are, therefore, the ratios giving zero mass-flow error. For any 
given profile and effective duct Mach nxmiber, the -Mach number ratios 
were higher for the area-weighting and mass-flow-weighting methods than 
for the mass-derived method, which explains the excessive mass flows 
computed by these methods. 

In the case of the momentum confutations, none of the averaging 
methods gives the base value of the Mach number ratio that is required 
to exactly compensate for the static -pres sure error and reduce the mom- 
entum error to zero. It is possible to determine the necessary Mach 
number ratio, however, by equating equations ( 4 ) and (lO). With the 
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separation profile at an effective duct Mach number of 0.71 (K = l)^ for 
example^ the ratio of measured to effective static pressiire is 0.753 for 
all averaging methods (fig. 2). In order to compute the correct thrust^ 
the calculated average Mach number should be 1.34 times the effective 
value. From figure 4^ the actual Mach number ratio is less than this 
value for the mass -derived method^ and more for the other methods. 

It is shown in the section ANALYSIS that the mass-flow errors in- 
oj dicated in figure 5 can be eliminated by redefining the Mach mmiber and 

cn static pressure of the uniform flow. The Mach numbers required to ac- 

conplish this are presented in figure 7. 

The Mach number errors for both the mass-flow-weighting and area- 
weighting methods are greatly reduced as conpared with the original er- 
rors shown in figure 4, In the case of the area-weighting method, the 
calculated Mach numbers are still larger than the effective values. 
However, the calculated Mach numbers are now lower than the effective 
values for the mass-flow-weighting method. These results are consistent 
with the calculated-to-effective total-pressure ratios shown in figure 
3. It can be deduced from equation (l4) that the Mach number ratios of 
figure 7 will be Inversely proportional to these total -pressure ratios. 

The values of static pressure required to satisfy the mass flow are 
compared with the effective values determined from the exact weighting 
procedure in figure 8. The errors in total-pressure calculation (fig. 

3) and Mach number calculation (fig. 7) tend to conp)ensate (eq. (15)), 
so that the static -pres sure error is greatly reduced as compared with 
the measured pressure shown in figure 2. The calculated static pres- 
sures for the area-weighting method were less than the effective value. 
For the mass-flow-weighting method, the calculated pressures exceeded 
the effective value. These trends arise from the predominant effect of 
total pressure, as compared with Mach number, in the static -pressure 
calculation. 

The momentums calculated with the static pressures and Mach numbers 
that satisfied the integrated mass flow are shown in figure 9. In gen- 
eral, these values are less in error than the values computed from the 
measured static pressure (fig. 6). An exception occurred with the area 
weighting of the separation profile. In this case, the calculated mo- 
mentum obtained with the assumption of measured static pressure was 
slightly greater than the integrated value, whereas that obtained for 
conditions satisfying the mass flow was less than the integrated value. 

The seriousness of the errors introduced by the various weighting 
methods depends on the use to which the averaged flow quantities are 
applied. The simple determination of diffuser total-pressure recovery, 
for exanple, is only subject to the errors indicated in figure 3. In 
the usual range of duct Mach number for which such data are evaluated 
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(Mach n-umbers less than 0.4), the errors associated with any of the 
weighting methods are small for the profiles examined. When the averaged 
quantities are to he utilized in broader applications, however, the er- 
rors arising from the various weighting methods may become more 
significant. 

Diffuser characteristics . - The diffuser pressure -recovery - air- 
flow characteristics that would be predicted by the varioias weighting 
methods for the separation profile are indicated in fig\xre 10. In the 
calculation of this figure, the average static pressure was assumed to 
correspond to the measured value. It was further assxmied that the ef- 
fective total-pressure recovery was 0.90 at an effective duct Mach num- 
ber of 0.3, corresponding to critical flow, and was constant in the sub- 
critical flow region of the inlet. ^ 

Inasmuch as the mass-derived method of averaging has no mass-flow 
error, the only difference between the diffuser characteristic pre- 
dicted by this method and the mass-momentum characteristic occurs in 
the level of the critical and subcritical pressure recoveries. The mass- 
flow errors Introduced by the mass -flow- weighting and area-weighting 
methods combine with the total -pressure errors associated with these 
averaging methods to cause marked shifts in the predicted diffuser char- 
acteristic as compared with the mass-momentum characteristic. In the 
supercritical flow region the corrected air flows predicted by the ap- 
proximate averaging methods at a given level of pressure recovery are 
in error in the same proportion as the mass -flow error indicated in 
figure 4. Conversely, at a given value of corrected air flow, large 
apparent differences in total-pressure recovery result with the various 
averaging methods. The choice of averaging method wo\ild thus have a 
large influence on the selection of inlet size to match a desired en- 
gine air -flow rate or on the prediction of the operating pressxcre- 
recovery level of an engine-inlet combination. 

The shift in apparent diffuser characteristic illustrated by fig- 
ure 10 would be less marked with the other profiles considered in this 
analysis, inasmuch as the total -pressure and mass-flow errors are 
smaller than for the separated profile. For weighting methods in which 
the mass flow is satisfied, the error in diffuser characteristic would 
be confined to the subcritical pres sure -recovery level, regeurdless of 
the profile. 


^Subcritical flow is defined as the regime where the absolute mass 
flow varies with changes in discharge pressure. When mass flow is in- 
dependent of back-pressure changes, the inlet flow is said to be super- 
critical. This is the hyperbolic region of the curves in figure 10. 

The intersection of these two flow regimes is termed the critical flow 
condition. 
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Diffuser drag calculations . - Figure 11 schematically illustrates 
the type of research model installation frequently used to evaluate 
combined internal and external flow problems of engine-inlet installa- 
tions . Although shown as a nose or nacelle inlet, the same type of in- 
stallation and support system can be used to study fuselage inlets. 

The model is supported from a sting by a balance which measvires the sum 
of the thrust and drag forces exerted on the model. Internal air -flow 
conditions are regulated by a plug in the discharge duct which is s\ip- 
ported from the sting. (This plug is generally remotely actmted to 
vary the air -flow conditions.) Internal flow conditions are evaluated 
by measurements at a flow measuring station in a region corresponding to 
the compressor inlet in the model prototype. 

Since the duct is cylindrical downstream of the flow measuring sta- 
tion, the only axial force on this section is a small viscous shear 
force which is generally neglected. The momentum evaluated at the meas- 
\aring station can therefore be used to determine the thrust force on the 
model. By subtracting the thrust force from the balance force, the ex- 
ternal drag of the model can be determined. 

It is shown in equation (B6) of appendix B that errors in momentum 
or mass -flow calculation at the measuring station cause errors in a 
drag-coefficient parameter according to the relation 



The magnitudes of these errors for a free-stream Mach nxamber of 2.0 are 
indicated in figure 12 for the various profiles and for weighting meth- 
ods in which the measured static pressure is satisfied. The sign con- 
vention is such that positives errors correspond to calculated drag coef- 
ficients that are less than the correct values. 

With each of the weighting methods, the error increased with in- 
creases in the duct Mach number in accordance with the increasing errors 
in mass -flow and momentum shown in figiires 5 and 6. In general, the 
mass-derived method, in which the mass flow as well as the measured 
static pressure is satisfied, gave the lowest drag errors. 

The importance of the errors indicated in figure 12 depends \ 5 )on 
the relative importance of the induction system to the over -all model. 
If, for example, the model represented by these error curves is a 
nacelle in which the duct area is 90 percent of the frontal area and 
the pressure recovery is 0.8, then the absolute error in drag coeffi- 
cient based on the frontal area would be 72 percent of the indicated 
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parameter error. At Mach number 2.0, the nacelle drag coefficient may- 
be on the order of 0.1 to 0.15 for an effective duct Mach ntmiber of 0.2. 
The indicated errors may thus become a large fraction of the desired 
value. If, on the other hand, the error ctarves of figure 12 apply to 
an inlet mo-unted on a fxiselage in which the duct area is a smaller frac- 
tion of the fuselage frontal area, the relative importance of the indi- 
cated errors is greatly reduced. 

The magnitude of the drag-coefficient errors due to errors in the 
weighting method is greatly reduced for weighting methods in which the 
integrated mass flow is satisfied, as shown in figure 13. The error 
curves for the mass-derived method are reproduced from figure 12. Both 
the mass -flow -weighting method and the area-weighting method produce 
less error than the mass-derived method with this criterion. As com- 
pared with the method in which the meas-ured static pressure was used 
(fig. 12), the errors introduced by the mass -flow -weighting method are 
reduced about 90 percent. For the area-weighting method, the errors 
with the mass flow satisfied are only on the order of one-fifth the er- 
rors when the measured static pressure was used. The sign of the errors 
obtained from the area weighting method is generally reversed between 
figures 12 and 13. This corresponds to the shift in value of the cal- 
culated momentxim relative to the true momentum shown between figures 6 
and 9. Except for the separation profile, the lowest drag-coefficient 
errors are obtained with the area-weighting method when the mass flow 
is satisfied. 

Figiores 12 and 13 illustrate possible drag-coefficient errors at a 
free-stream Mach ntanber of 2.0. The effect of free-stream Mach nvmiber 
is illustrated in figure 14. The profiles evaluated in this fig\ire all 
have a maximum duct Mach nvmiber of 0.4, which corresponds to an effec- 
tive Mach nvmiber of about 0.35 in each case. These calculations are 
for the weighting methods in which the uniform-flow static pressvire is 
assumed equal to the measvired value; hence, the mass-flow errors indi- 
cated in figure 5 are included. Similar trends would be observed for 
the weighting methods in which the integrated mass flow was satisfied. 
The magnitude of the drag -parameter errors would be decreased in the 
latter case, however. 

The increasing error in drag parameter with increasing svpersonic 
Mach number does not necessarily imply an increase in the absolute 
drag-coefficient error of the same proportion. The total -pres sure - 
recovery term in the denominator of the drag parameter will generally 
decrease with increasing Mach nvmiber. This will compensate in part for 
the increase in parameter error. For such cases, the anticipated error 
in drag coefficient may remain relatively constant throu^out the svper- 
sonic Mach nvmiber range. If, on the other hand, highly efficient in- 
lets are being considered at high Mach nvmibers, the drag-coefficient 
error will increase for a given level of flow distortion as conpared 
with the errors resulting at lower Mach nvmibers. 
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Inlet pressure recoveries may be expected to remain at a generally 
high level throughout the subsonic Mach number range. It would be antic- 
ipated that the potential error in drag coefficient would therefore in- 
crease as Mach number is reduced unless there was a concomitant improve- 
ment in the duct profile. 

Variable -area^^duct calculations . - In many duct flow applications, 
uniform-flow properties are calculated at a flow measuring station by 
one of the wei^ting methods. One -dimensional flow equations are then 
used to compute flow properties at other stations in the duct by the 
assumption of appropriate total -pres sure losses. These resultant prop- 
erties are affected by the errors previously demonstrated to be associ- 
ated with the various weighting methods. Additional errors are intro- 
duced if there are area changes in the duct. 

The nature of the errors introduced in variable-area-duct calcula- 
tions can be illustrated by the flow shown in figure 15. It has been 
assumed in this flow that a uniform static pressure exists at each sta- 
tion and that each filament of the flow expands isentropically between 
the two stations. 

Each filament diffuses to a higher static pressure as the flow pas- 
sage area increases. The static-pressure rise is constant across all 
filaments; consequently, the filaments having low initial velocity under- 
go a greater deceleration than those with high velocity. The expansion 
rate varies as a result, and the low-velocity filaments occupy a larger 
fraction of the final duct area than of the initial duct area. 

As previously shown, the mass, momentimn, and energy of the nonuni- 
form flow at each station can be duplicated by a uniform flow determined 
by the mass -momentum method. The resultant average total pressure at 
each station includes the mixing losses that would be incurred if the 
nonuniform flow were allowed to mix fully in a constant-area section. 

The magnitude of the mixing losses depends on the velocity differences 
between fluid filaments in the nonuniform flow. These differences are 
greater after diffusion than in the initial flow. Thus, the uniform 
flow satisfying the mass, momentum, and energy of the real flow must 
undergo an apparent total -pressure loss in the diffusion process, even 
though the real flow expands isentropically. A final flow calculated 
from the average initial flow by isentropic one -dimensional equations 
will therefore be in error. 

The magnitude of the errors introduced through the assumption of 
one -dimensional average flow properties is illustrated in figure 16. 

For this example, the initial profile was assumed linear with K = 1.0 
(eq. (l9)). The final profiles and duct areas were analytically deter- 
mined for a range of static -pres sure ratio for assumed Isentropic ex- 
pansion of the nonuniform flow by the method outlined in appendix C. 
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Average flow properties were determined at each station hy the conven- 
tional weighting methods as well as hy the mass-momentum method. The 
figure presents the ratio between the average weighted properties at 
each station and those calculated by applying isentropic one -dimensional 
relations to the initial weighted flow. 

As previously indicated, there is an effective loss in total pres- 
sure in the expansion process when evaluated by the mass-momentum method. 
Similar losses are calculated by the area-weighting and mass-derived 
methods. In addition to the loss in total pressure, the average values 
of Mach number and the calculated momentum and mass flow are lower at 
each station in the duct than would be predicted by the one -dimensional 
calculation. 

If the average flow properties at each duct station are determined 
by the mass-flow-weighting method, the one -dlmens ional equations may be 
applied without error. With this weighting method, each filament of the 
nonuniform flow exerts a weight in the average total -pres sure determi- 
nation that is proportional to its increment of mass flow and total 
pressure. These quantities remain invariant in the expanded filament; 
consequently, the calculated average total pressure remains constant. 

The error shown in figure 16 for each weighting method is a rela- 
tive error for the given expansion ratio. It represents the difference 
between the value of the flow property as con 5 )uted from one -dime ns ional 
relations and the value determined from a weighting of the local flow. 
The previously discussed inherent error between the weighted flow prop- 
erties and the integrated flow properties must also be considered be- 
fore the absolute error associated with the application of one- 
dimensional relations to variable-area duct flows can be determined. 


CONCLUDING REMARKS 

It has been shown that conventional weighting methods used to ob- 
tain uniform flow representations of nonuniform duct flows can cause 
large errors in the calculated uniform-flow properties. These errors 
are predominantly associated with the conventional assumption that the 
measured static pressure can be used in conjunction with a weighted 
total pressure to define the uniform flow. 

An averaging method has been developed which yields uniform-flow 
properties that reproduce the mass, momentum, and total energy of the 
nonuniform flow without error for special cases in which the total tem- 
perature and static pressure are constant across the duct. The magni- 
tude of the errors introduced by conventional weighting procedures may 
often warrant the additional complications required to apply this method. 
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It has also been shown that nonuniform flows through variable -area 
duct passages result in changes in average flow properties that are not 
associated with the real thermodynamic flow path. Consequently, addi- 
tional errors are introduced into nonuniform duct flow calculations when 
one -dimensional equations are applied to the averaged flow at one station 
in order to predict the averaged qtiantities at another station. 

These findings indicate that care should be exercised in the selec- 
tion of a method of averaging nonuniform duct flows and that calcula- 
tions based upon the weighted flow should be interpreted with caution. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 13, 1954 
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APPEOT)IX A 
SYMBOLS 

The following symbols are used in this report : 

A flow area 

A^gj reference area for drag coefficient 

(O 

aa stagnation speed of sound 

Cj) drag coefficient, Drag/qQAj.gj. 

net thrust 

K . maximvim duct Mach number 

M Mach number 

m mass -flow rate 

P total presstire 

p static pressure 

q dynamic pressure, q = ■|pV^ = ^M^ 

R gas constant 

T absolute total temperat\ire 

V velocity 

W weight-flow rate 

X fractional distance from wall to duct centerline 

X ratio of specific heats, 1.4 for air 

p mass density 

5 total pressure, corrected to NACA standard sea-level conditions, 

P/2116 
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9 total temperatvire , corrected to MCA standard sea-level conditions, 
T/519 

<P momenttun, <p = mV + Ap = pA(l + yM^) 


Subscripts : 


c calctilated 

OJ 

^ e effective 


i initial station in an expanding duct 


is isentropic 


0 free stream 
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APPENDIX B 


CALCULATION OF DRAG ERRCKS IN DUCTED-BODY INVESTIGATIONS 

The net internal force acting on a ducted body is the difference 
between the outlet and free-stream momentum. If the model is similar 
to that shown in figxare 11, in which the duct is cylindrical downstream 
of the force measuring station, the only axial force on this section 
will be a small viscous shear force. This shear force is generally 
neglected, and the momentum evalmted at the measuring station is as- 
simied equal to the outlet momentum. The net internal force therefore 
becomes 


Fq = - ^0 (bi) 

The absolute error in net thrust arising from errors in the deter- 
mination of the momentum and mass flow in the duct becomes 



where <p and m are the integrated values of momentum and mass flow, 
respectively, and and m^, are calculated values based upon in- 

exact averaging methods. 

The absolute qmntities in the terms on the right side of equation 
(B2) can be reduced to functions of the equivalent duct Mach number by 
introducing the measured duct static pressvire and the total temperature, 
which gives 


pA 



(B3) 


Since the balance measures the sum of the thrust and drag forces 
on the model, the error in calculated external drag will be numerically 
equal to the error in calculated thriist from equation (B3). The re- 
sultant error in drag coefficient based on any arbitrary reference area 
is 


, _ '^n _ 2 A p ^n 

~ Wef tM^ \ef PO PA 


(B4) 
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By using equations (B3) and (B4) and ttie relation 

_p p_ £e ^ ^ 

Pq ■ Pe Pe PO ^0 


(B5) 


the following drag -coefficient -err or parameter can he determined, which 
is a function of free-stream and measuring-station flow conditions only: 


ACt 


PO Aref 


2 Pq p Pe 
^0 ^e ^e 


JL 

pA 


RT 


(B6) 
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APPENDIX C 


DETERMINATION OF NONUNIFCKM-FLOW PROFILE AFTER ISENTRCPIC DIFFUSION 
The continuity equation may be written in differential form as 


dm 





dA 


(Cl) 


Thus, 


»l(l ^ <^1 


(C2) 


where the subscript i refers to the initial duct station before 
diffusion. 


For isentropic flow, 


.P. = E ^ 

Pi p Pi 


-1 



Combining equations (C2) and (C3) gives 


(C3) 


dA = 


Mi dA^ 



r-1 


1 -l^ 

i. 


ni‘ 

r 


71^ 


(C4) 


For two-dimensional flow, dA^/A^ = dXj^, The required flow area 
after diffusion therefore becomes 


9622 


NACA TN 3400 


23 



The Mach number of any filament after diffusion is, from equation 

(03), 



The coordinate x of the filament after diffusion is 



In the exan^jle considered herein, the initial profile was assumed 
to follow the linear equation = 0.2Xj^ + 0.8. From equation (C5) the 

required flow area after diffusion becomes 



The flow coordinate for a given filament of the flow becomes 
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r r-i] r y-r 

1 -/ ^ - A 0. 64 + 1 - /-£-\ ^ 

r-1 [p±J ] r-i [pj 


1 + -^ 1 -(^) -/Jo. 64 + -^ 1 -f— 

r-1 [p±J i r-1 IPi 


Upon substitution for from equation (C3), equation (C9) may- 

be simplified to 


t ^ -[fl 


^ 1 j£-\ >■ 




(h) 



r-1 

2 

/p\ 

0.64 + -Ar 

1 - 

r-1 

\Piy J 


(CIO) 


The profile after diffiosion is therefore also linear, and the 
■weighting eqTiations may be solved directly for the uniform-flow 
properties . 
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Figure 1. - Typical profil 
duct Mach number, 0.2. 
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(c) Area-weighting method. 


Figure 3 
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(a) Mass-flow- weighting method. 
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(b) Area-weighting method. 


Figure 5. - Mass flow calculated by two weighting methods 
Static pressure equal to measured value. 
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(a) Mass-derived method 
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(b) Mass-flow- weighting method 


Effective duct Mach number, Mg 

(c) Area-weighting method. 

Figure 6. - Momentum calculated by three weighting metho. 
Static pressure equal to measured value. 
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(a) Mass-flow- weighting method. 
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Effective duct Mach number, M 
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(b) Area-weighting method. 

Figure 7 . - Duct Mach numbers calculated by two weighting 
methods: integrated mass flows satisfied. 
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Figure 15. - Typical profile distortion in an isentropic diffusion process. 
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Figure IG. - Comparison of locally weighted and isentropically calculated flow 
properties in an expanding duct. 
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